JOURNAL OF RESEARCH of the National Bureau of Standards 
Vol. 89, No. 3, May-June 1984 



A Head-Space Method for Measuring 

Activity Coefficients, Partition Coefficients, and Solubilities 

of Hydrocarbons in Saline Solutions 

Stanley P. Wasik, Frederick P. Schwarz, Yadu B. Tewari, and Michele M. Miller 

National Bureau of Standards, Gaithersburg, MD 20899 

and 

J. H. Purnell 

University College of Swansea, Swansea, Wales, UK 



Accepted: February 2, 1984 

An apparatus is described which measures the equilibrium distribution of a hydrocarbon between a gas and 
aqueous phase. Soluble hydrocarbons are extracted from an aqueous salt solution by very small bubbles of 
hydrogen generated electrolytically from a gold electrode located at the bottom of a cylindrical cell. The 
partition coefficient is determined from the volume of the aqueous solution and the solute concentration in the 
head-space after a measured volume of hydrogen has bubbled through the cell. The concentration of the solute 
in the head-space is measured by gas chromatography. The observed distribution is supplemented by vapor 
pressure and molar volume data and can be used to calculate the solubility and the activity coefficient of the 
solute in the aqueous phase. The partition coefficient, activity coefficient, and solubility for 18 alkylbenzenes in 
aqueous 0.5 M NaCl at 25 °C were measured by this method. 
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TntrnrinHinn search, it is useful for understanding how hydrocarbons 

migrate and accumulate to form oil deposits [3]. In biol- 

A knowledge of the equilibrium properties of aqueous ogy, a knowledge of ho w hydrocarbons behave in aque- 

hydrocarbon solutions is valuable in several fields. In ous solutions is important for understanding the effects 

water pollution control, such information is helpful in of hydration on the configuration of biopolymers [4], In 

devising abatement processes [l] 1 , in modeling natural chemistry, experimental data on these systems are 

water systems [2], in designing toxicity experiments, and needed for testing models of water and aqueous solu- 

in developing analytical methods. In petroleum re- tions [5]. 

These aqueous solutions can be characterized by de- 
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saturation vapor pressure. Another thermodynamic 

1 Figures in brackets indicate literature references at the end of property of interest is the solute activity coefficient 

this paper. based on volume fraction, y^, which may be calculated 
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from the value of K, the solute molar volume, and the 
solute saturation vapor pressure. This property is of par- 
ticular importance in developing correlations based on 
some additive property of the solute, such as molar vol- 
ume or carbon number, since y$ depends only weakly on 
temperature and varies significantly less with solute 
structure and configuration than does the corresponding 
correlation with solubilities or with the mole fraction- 
defined activity coefficient, y x . 

We have previously described a method for mea- 
suring K using a head-space technique [6, 7, 8] which is 
best suited for measuring K in the range 1 — 10. We now 
report a method without limitations on the value of K. 
In this method, the soluble hydrocarbons are extracted 
from aqueous salt solutions by very small bubbles of 
hydrogen generated electrolytically from a gold elec- 
trode located at the bottom of a cylindrical cell. The 
partition coefficient is determined from the volume of 
aqueous solution and the solute concentration in the 
head-space after a measured volume of hydrogen has 
bubbled through the cell. The concentration of the hy- 
drocarbon in the head-space is analytically determined 
by conventional gas chromatography. The partition co- 
efficients, solubilities, and activity coefficients for 18 
alkylbenzenes in aqueous 0.5 M NaCI at 25 °C measured 
by this method are reported here. 



Experimental 

Figure 1 shows the extraction cell. The compartments 
A and B have internal volumes of 50 and 40 ml, re- 
spectively. The two compartments are connected by 
means of a 14 mm o.d. glass tube with a 10 mm o.d. 
coarse porosity glass frit on one end. A silicic acid plug 
is precipitated on the compartment B side of the frit to 
prevent the flow of liquid from one compartment to 
another and to provide low resistance for the flow of 
electric current. The gold electrode (1 mm o.d. wire) in 
compartment A is positioned by means of a septum held 
in a l/4x 1/4 in union fitting. This electrode has a spiral 
configuration to provide a large surface area. The elec- 
trode in compartment B is a straight segment of 1 mm 
o.d. platinum wire. 

The extraction cell is immersed in a bath controlled to 
±0.02 °C. The thermometer used in this work was 
calibrated with an accuracy of ±0.01 °C by the National 
Bureau of Standards. The head-space is sampled by a gas 
injection valve thermostated at 150 °C. The 1/8 in o.d. 
stainless steel tube connecting the extraction cell with 
the valve is heated to 1 50 "C by a heating tape. 

The chromatographic column is a 15 mX 0.5 mm i.d. 
SCOT column prepared with finely ground di- 
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Figure 1 — Extraction cell compartments A and B have internal 
volumes of 50 and 40 ml, respectively. 



atomaceous earth on a fused silica support and coated 
with a mixture of m-bis(m-phenoxy phenyl) benzene and 
Apiezon L. The effluent is monitored by a hydrogen 
flame detector with an electronic integrator measuring 
peak areas. 
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A constant current power supply with a stability of 
1/10,000 provides electrical current to the extraction 
cell. The volume of hydrogen plus water vapor, V, bub- 
bled through compartment A is calculated from the 
equation 



V={it/2F) {RT/{p A -p w y\ 



(1) 



where / is current; /, the time in seconds; F, the Far- 
aday constant; R , the gas constant; p\, the atmospheric 
pressure; and^v the vapor pressure of the aqueous solu- 
tion at temperature T(K). 

The aromatic hydrocarbons used in this study had a 
stated purity of 99.9 mole %. The water used in the 
preparation of solutions was doubly distilled over potas- 
sium permanganate. 



Procedure 

A mixture of several of the alkylbenzenes was pre- 
pared such that their partial pressures at 25 °C were 
approximately the same. Five ml of the mixture was 
gently stirred with 100 ml of aqueous 0.5 M NaCl for 2 
h and then allowed to stand overnight. Compartment A 
was filled with the aqueous solution leaving only a small 
volume of head-space (< 1.0 ml) above the solution. The 
volume of the solution was determined from the weight 
of the cell before and after filling and from the density of 
the solution. Compartment B was filled to the same level 
as compartment A with 5% v/v H 2 S0 4 solution. The 
extraction cell was then immersed in the constant tem- 
perature bath for 1 h before electric current (0.2—0.4 
amp) was passed through the cell. A small volume of the 
head-space (0.3 mL) was injected into the gas chro- 
matograph via the gas sampling valve for analysis. The 
alkylbenzenes were chosen so that there was a baseline 
separation for each peak. The peak areas were recorded 
along with the time of sample injection. This procedure 
was repeated every 10 min until the peak areas were 
approximately 1/100 their original value. 

In order to measure the peak area corresponding to 
the solute saturation vapor pressure three columns were 
constructed of stainless steel tubing, 2.5 mm i.d. and 0.5 
m in length. The solid support was Chromosorb W 
(60—80 mesh) with toluene, ethylbenzene, and n- 
propylbenzene as the stationary phase, respectively. 
Each column was immersed in the constant temperature 
bath, and helium saturated with water vapor was al- 
lowed to flow through it before being sampled for anal- 
ysis. 



Thermodynamic Background 

Consider a known volume, V L , of aqueous salt solu- 
tion of hydrocarbons in a cylindrical vessel containing a 
gold electrode at the bottom. Connected to this vessel 

by means of a glass fritted disc is another vessel contain- 
ing an aqueous solution of H 2 S0 4 . When electric current 
with the appropriated polarity passes through the cell, 
hydrogen in the form of very small bubbles is evolved at 
the electrode in the vessel containing the aqueous solu- 
tion, and oxygen is evolved at the other electrode. Wa- 
ter vapor and the dissolved hydrocarbons in the aqueous 
solution are equilibrated with the hydrogen bubbles. 
The hydrocarbon concentrations in the head-space are 
measured by means of gas chromatography. We may 
assume that the concentration of a particular hydro- 
carbon in the head-space is proportional to its concen- 
tration in the aqueous solution. One then has the re- 
lationship 



CJy)=K l C u (y) 



(2) 



where C wi (v) and C hi (y) are the concentrations of the /' h 
hydrocarbon in the aqueous solution and in the head- 
space, respectively, after a volume V of hydrogen plus 
water vapor has bubbled through the cell and K f is the 
proportionality factor. Under equilibrium conditions K t 
would be the partition coefficient. Conservation of mass 
requires that the ratio of C hi (V), corresponding to the 
passage of volume V, to its value, C hi (0), at zero hydro- 
gen volume is given by the expession 



In [C„( V )/C„(0)] = - V/(V L K,) 



(3) 



The peak area for the /'* hydrocarbon, A t {V) may be 
expressed as 



A,(V) = cr l C u (V)V, 



(4) 



where cr, is the instrument sensitivity constant for the I th 
hydrocarbon and V s is the volume of the gas sampling 
valve sample loop. Substituting for C hi (V) in eq (3), 



inN,(r)/4(0)]=-r/r L j:,. 



(5) 



A plot of In At(V) versus V (or time) should thus be 
linear if K t is constant, and the slope will give the value 
ofK,. 

The attainment of equilibrium of the dissolved hydro- 
carbons in the aqueous solution with the hydrogen bub- 
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bles depends on the size of the bubbles and the contact 
time. Hydrogen bubbles generated electrolytically are 
less than 0. 1 mm in diameter. Under these conditions the 
bubbles are small enough and the contact time with the 
solution long enough that equilibrium would be ex- 
pected before the bubbles reached the surface, and the 
measured K t should be the partition coefficient. 

During the course of the measurements water is con- 
tinually being removed from the generating cell as satur- 
ated water vapor. For measurements near room tem- 
perature the amount of water removed is very small 
( < 0. 1 % ); hence, no corrections are made for this effect. 



Activity Coefficient 

The partition coefficient K t may be defined as 



K t = 



C wi (y 
C M (y) 



)~\Vl) Pi 



(6) 



where n, is the number of moles of the /* hydrocarbon, 
P t is the partial pressure of hydrocarbon at temperature 
T, R is the gas constant, and V L is the total volume of 
aqueous solution. The partial pressure (P,) may be ex- 
pressed in terms of volume fraction-based activity coef- 
ficient iy^) as 



P,. ^y^PI 



(7) 



where P° is the saturated vapor pressure of the solute 
and 4>i is the solute volume fraction combining eqs (6) 
and (7) we get 



Ki = 



RT 

VvViP", 



(8) 



where V, is the molar volume of the z ,h hydrocarbon at 
temperature T. 



Solubility 

In order to calculate the hydrocarbon aqueous solu- 
bility from K values, it is essential that the hydrocarbon 
concentration in the aqueous phase be proportional to its 
concentration in the gas phase, P t /RT, up to the solute 
saturation vapor pressure, pi, i.e., a two-phase system, 
liquid solute and aqueous solution saturated with solute. 
Under these conditions 



where CI is the solubility of the solute temperature T. 
When liquid hydrocarbon is present at equilibrium with 
an aqueous phase, it will dissolve some water. The 
amount of water dissolved is so small, however (the 
solubility of water in all hydrocarbons at 25 "C is well 
under x =0.01) that no correction need be made for its 
effect on the hydrocarbon's vapor pressure. 



Results and Discussion 

Ln A(V) versus V plots (eq (5)) for some al- 
kylbenzenes in which the initial aqueous solution was 
saturated with the solute indicate that the propor- 
tionality between solute concentration in the gas phase 
and aqueous phase holds up to the solute saturation 
concentration in the aqueous phase (solubility) and that 
eq (9) may be used to calculate solubilities for these 
compounds. Values for ln A t (0) (time equal to zero) 
were obtained from the saturation vapor pressure at 
25 "C using pure solute. 

The partition coefficients were calculated from the 
slopes of the \nAi{V) versus V plots for the al- 
kylbenzenes. These values are given in table 1 along 
with values of the solute activity coefficients calculated 
via eq (8). For comparison, values of the solute activity 
coefficients calculated from solubility data measured us- 
ing the generator column method [10] are given in the 
same table. The good agreement between the two sets of 
data justifies the assumption made in deriving eq (7), i.e., 
■y^l for an apolar solute in equilibrium with an aque- 
ous phase. 

The head-space method proposed in this paper may 
be used to measure K for any type of compound. The 
largest error in the measurements comes from mea- 
suring peak areas. Since one needs to know the peak 
area at time t , as well as the peak area corresponding to 
the saturation vapor pressure, the average error in the 
partition coefficients is approximately 1.5%. 

The method has a disadvantage in that K can be mea- 
sured only in salt solutions since some electrolyte is 
required in the solution to conduct the electrical cur- 
rent. In order to determine K in pure water, K must be 
measured in salt solutions of differing concentration, C s , 
and then extrapolated to zero salt concentration using 
Setchenow's expression 



log(tfVAT™)=^Q 



(10) 



K=ClRT/p° 



(9) 



where K w h the solute partition coefficient in water, K m 
is the solute partition coefficient in salt water of concen- 
tration C s , and k, is the salting-out coefficient. 
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Table 1. Partition coefficients, solubilities and activity coefficients for alkylbenzenes in aqueous 0.5 M NaCl at 25.0 °C. 





Pi 


V, 


K.a/w 


cyxio' 


logy* 


Solutes 


(mm) 


(mL/mole) 




(M) 




Benzene 


95.18 


89.4 


3.47±0.03 


17.76 


2.799 


Toluene 


28.45 


106.9 


2.74±0.02 


4.19 


3.349(3.292)' 


Ethylbenzene 


9.51 


123.1 


2.19±0.02 


1.12 


3.861(3.809) 


1 ,2-Dimethylbenzene 


6.62 


121.2 


3.52 ±0.02 


1.25 


3.820(3.735) 


1 ,4-Dimethylbenzene 


8.84 


123.9 


2.24 ±0.02 


1.06 


3.882 


n-Propylbenzene 


3.43 


140.1 


1.53±0.01 


0.282 


4.403(4.419) 


Isopropylbenzene 


4.64 


140.2 


1.20±0.001 


0.299 


4.378(4.260) 


1,3,5-Trimethylbenzene 


2.61 


139.6 


2.01 ±0.01 


0.282 


4.405 


1,2,4-Trimethyibenzene 


2.096 


137.9 


2.78±0.02 


0.313 


4.365 


2-Ethyl-l-Methylbenzene 


2.521 


137.1 


2.66±0.02 


0.361 


4.305 


3-Ethyl- 1-Methylbenzene 


2.999 


139.7 


1.69±0.02 


0.273 


4.419 


n-Butylbenzene 


1.011 


156.8 


1.18±0.01 


0.0642 


4.997(4.952) 


Isobutylbenzene 


1.929 


158.1 


0.677+0.01 


0.0702 


4.955 


sec-Butylbenzene 


1.873 


156.4 


0.870+0.01 


0.0876 


4.863 


t-Butylbenzene 


2.208 


155.6 


1.13±0.01 


0.134 


4.681 


1 ,3-Diethylbenzene 


1.134 


156.1 


1.22±0.01 


0.0744 


4.935 


1 ,2-Diethylbenzene 


1.043 


153.2 


1.81 ±0.02 


0.1015 


4.794 


1 ,4-Diethylbenzene 


1.051 


156.4 


1.26±0.01 


0.0712 


4.953 



The values in brackets are calculated from experimental solubility data using the equation y t =(C7 F<) -1 where P,is the molar volume of 
the solute [10]. 
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